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Under certain conditions, sound waves in a fluid may be governed by a Klein-Gordon equation 
on an 'effective spacetime' determined by the background flow properties. Here we consider the 
draining bathtub: a circulating, draining flow whose effective spacetime shares key features with the 
rotating black hole (Kerr) spacetime. We present a complete investigation of the role of quasinormal 
(QN) mode and Regge pole (RP) resonances of this system. First, we simulate a perturbation in the 
time domain by applying a flnite-difference method, to demonstrate the ubiquity of 'QN ringing'. 
Next, we solve the wave equation in the frequency domain with the continued- fraction method, to 
compute QN and RP spectra numerically. We then explore the geometric link between (prograde 
and retrograde) null geodesic orbits on the spacetime, and the properties of the QN/RP spectra. 
We develop a 'geodesic expansion' method which leads to asymptotic expressions (in inverse powers 
of mode number m) for the spectra, the radial functions and the residues. Next, the role of the 
Regge poles in scattering and absorption processes is revealed through application of the complex 
angular momentum method. We elucidate the link between the Regge poles and oscillations in the 
absorption cross section. Finally, we show that Regge poles provide a neat explanation for 'orbiting' 
oscillations seen in the scattering cross section. 

PACS numbers: 04.70.-s, 04.30. Nk, 43.20.-l-g, 47.35.Rs, 11.55.Jy 



I. INTRODUCTION 

Black holes - 'trapped' regions of spacetime - are a 
key element of Einstein's theory of general relativity. Al- 
though perhaps once viewed as mathematical curiosi- 
ties, astronomers have now compiled a range of com- 
pelling evidence for their existence. Black holes are a 
key ingredient in modern theories of galaxy formation, 
quasars, accretion disks, gamma ray bursts and super- 
novae. Yet, even if black holes were nothing more than 
a theorist's 'thought experiment', they would still have 
provoked the development of theoretical physics. In the 
1970s, Hawking and others showed that quantum field 
theory in curved spacetime implies that black holes are 
not completely black: they must radiate thermally, with 
a negative heat capacity. Furthermore, black holes seem 
to have a well-defined entropy, that scales with the hori- 
zon area [T]. This realisation has inspired myriad (and 
ongoing) attempts to consistently combine relativity and 
field theory in the strong field regime. 

There is more to black holes than Hawking radiation. 
A key property of a black hole is that it bends and traps 
light. Light rays may orbit a black hole in the vicin- 
ity of a 'photon sphere' which lies somewhat outside the 
horizon (at r = 3r/i/2 for the Schwarzschild BH, where 
is the horizon radius). The existence of an unstable 
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photon orbit gives rise to various interesting effects, such 
as the strong-field lensing of light from a distant source 
passing close to a black hole. In particular, the photon 
orbit is intimately linked to characteristic 'damped res- 
onances' that appear when waves interact with a black 
hole. Mathematically, damped resonances are manifest 
as poles in the scattering matrix S. The poles occur at 
complex frequencies and at complex angular momenta, 
and the corresponding modes are known as quasinormal 
(QN) and Regge pole (RP) modes, respectively [2H1]. In 
this paper, we investigate QN and RP modes and their 
physical consequences, in the setting of a simple rotating 
'black hole analogue' [7]. 

It seems unlikely that we will ever study black holes 
directly in the laboratory. Yet, as Unruh [5] noted three 
decades ago, we may study analogues: artificial systems 
(in various media) which exhibit some key kinematic fea- 
tures of black holes . For example, sound waves in fluid 
flows which are inviscid, irrotational, and barotropic are 
governed by the same wave equation as a scalar field in 
a curved space-time, namely [7] 

= -^9^ (VWla'^'d,^) = 0, (1) 

where here g^j^ is the effective metric (with inverse g^^ 
and determinant g). Note that in this context, g^j^y de- 
pends algebraically on the local properties of the fluid 
flow, and it does not (in general) represent a solution of 
the Einstein equations. Nevertheless, it is an intriguing 
prospect that, by studying sound waves on a background 
flow, one may understand better the propagation of fields 
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on a curved spacetime. In recent years, a wide range of 
black hole analogues in various media have been pro- 
posed and, indeed, studied, in the laboratory [TDHH]- A 
surge of recent experimental activity appears to be bear- 
ing fruit, as evidenced by a recent claim of experimental 
observation of correlations related to Hawking radiation 
in a wave-tank [15j. In this experiment, instead of a 
black- hole analogue, a white-hole analogue is used. A 
further example of a simple 'white-hole' analogue in flu- 
ids is the so-called circular hydraulic jump [16]. 

One of the simplest analogue models is the so-called 
draining bathtub (or draining vortex), described in [7]: a 
two-dimensional circulating flow with a sink at the ori- 
gin. In 2002, Schiitzhold and Unruh [17] described a 
possible experimental realization of the 'bathtub' idea, 
wherein gravity waves propagate in a flowing fluid in 
a shallow basin of varying height h{r). Non-dispersive 
long-wavelength perturbations are governed by an effec- 
tive geometry with line element ds^ = g^^dx^dx'^ , where 



ds' = -c'dP 



DdiV ( ~ Cdi'' ^ 

dr H + rd(j) 

r J \ r 



(2) 



Here the constants of circulation (C) and draining (D) 
[l8] relate to the background flow velocity vq of the fluid, 
namely 



Vq = —Dr/r + C(j)/r. 



(3) 



We assume that D > 0, so that the system acts as a 
black, rather than white, hole. In the model described in 
Ref. , the speed of the perturbation in the fluid c is set 
by = ttghao, where ag is the acceleration due to grav- 
ity and /loo is the height of the fluid far from the centre. 
Note that we assume we are within the linear dispersion 
regime, so that the perturbations propagate with a con- 
stant speed (i.e. c is independent of frequency). The ana- 
logue event horizon (where the inward flow rate exceeds 
c) lies at rh = -D/c, and the analogue ergosphere (where 
the flow becomes supersonic |vo| > c) has a boundary 
at Te = -y/C^ -I- D^/c [7]. Small perturbations 5w to the 
flow, V = Vq -l- (Sv, may be expressed in terms of a gra- 
dient of a potential 5-v — — V<I>, and the potential field 
<I> satisfies the Klein-Gordon equation ([T]) with effective 
metric ^ . Henceforth we set the speed of the perturba- 
tion equal to unity (c = 1). 

There are several key motivations for considering the 
draining bathtub (DBT) . Firstly, as described above, the 
DBT may perhaps be realized in the laboratory. Sec- 
ondly, the DBT provides a useful 'toy model' for the 
most astrophysically-relevant black hole, i.e. the Kerr so- 
lution. For example, both the DBT and the Kerr solution 
possess a horizon and an ergosphere, and can exhibit su- 
perradiance; but whereas the angular momentum of the 
(non-naked) Kerr black hole is constrained, J < Af^, 
the angular momentum of the DBT is (in principle) un- 
bounded. This follows as a consequence of their differing 
symmetry: the DBT is cylindrically symmetric, whereas 
the Kerr solution is axially symmetric. A third reason is 



simplicity: the DBT is arguably the simplest asymptot- 
ically flat rotating spacetime that can be envisaged (see 
also the cosmic string [IH])- It serves as a testing ground 
for developing calculation methods that can be extended 
to Kerr geometry. 

Given such motivations, it is no surprise to find that 
the spectrum of quasinormal modes of the DBT has al- 
ready received some attention [501 HI]- On the other 
hand, the Regge pole spectrum has not been considered. 
Given the close relationship between QN and RP modes, 
we believe that a comprehensive study of the resonances 
of the DBT is now justified. Here we can go many steps 
beyond existing work to show: (i) how quasinormal res- 
onances arise in time-domain simulations of a small per- 
turbation in the flow, (ii) how QN and RP resonances are 
closely related to the properties of the co- and counter- 
rotating null orbits, (iii) how complex angular momen- 
tum methods [52H23] may be applied to compute absorp- 
tion and scattering cross sections [551 US j (iv) how the 
fine-structure of the absorption cross section is related to 
the Regge pole spectrum [27l |28] , and (v) how 'orbiting' 
oscillations in the scattering cross section are also related 
to Regge pole spectrum. 

The remainder of this paper is structured as follows. In 
Sec. |TT] we study perturbations and null geodesies in the 
DBT effective spacetime. In Sec. |III| we evolve Gaussian 
initial data in the time domain and identify the quasi- 
normal mode ringing signal. We apply the continued- 
fraction method to obtain the frequency spectrum of the 
DBT QN modes. We extend the geodesic expansion 
method of Ref. [29l[30] and apply it to find the QN modes 
of the DBT. We validate the expansion method formulae 
by comparing with numerical results from the continued- 
fraction method. In Sec. IV we extend the expansion 



method to find approximations for the Regge poles, and 
again check against numerical data. In Sec.|V]we harness 
the power of the Complex Angular Momentum method 
to understand the key features of the DBT absorption 



and scattering cross sections. In Sec. V A we show that 
the absorption cross section CTabs('^) of the DBT [25] may 
be expressed in terms of the Regge poles |271 US] , and we 
derive a simple 'geometric' approximation for the cross 
section. In Sec. |VB| we show that 'orbiting' oscillations 
which arise in the scattering cross section are also re- 
lated to the Regge poles. We extend the geometric ex- 
pansion method (in Appendix A) to obtain an approxi- 
mation for the residues of the Regge poles, which we use 
to derive a simple geometric approximation for 'orbiting' 
at large scattering angles. We validate our approxima- 
tions against numerical data. We conclude with our final 
remarks in Sec. I VII 



II. THE DRAINING BATHTUB 

The draining bathtub was briefly described in the pre- 
vious section, and we address the reader to Ref. [7] for 
more details. In the lab-based coordinates {t,r,(j)}, the 
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effective geometry is described by the line element 
Following [20], it is convenient to introduce an alterna- 
tive coordinate system {t,r, (p} via dt — di ~ Ddr/{rf), 
dcj)^ d(j) — CDdr/{r^f), with (/)(r — > oo) — > and 

f{r) = l-Dyr\ (4) 

The line element ^ then takes the form 

ds^ = -f{r)df + f{r)-^dr^ + {rd(t) - Cdt/rf . (5) 

Henceforth we will work with these (non-lab based) co- 
ordinates exclusively. 

A. Perturbations 

Small perturbations to the background flow, Sv — 
— V$, are governed by the Klein-Gordon equation ([I]) 
with line element ([5|. Let us now decompose $ in az- 
imuthal modes, namely 



$(i,r. 



^ OO 



im(p 



(6) 



m— — OO 



Since $ is a real field, the symmetry relation — ip^m 
follows. Inserting Eq. ^ into Eq. ([T]) leads to the wave 
equation 



d iCm\ 9^ , s 



^^{t,r) = 0, (7) 



where 



4r4 



(8) 



and the tortoise coordinate is defined by dr^ = / ^dr, 
or, explicitly. 



D 



In 



r- D 



r + D 



(9) 



A perturbation of compact support in the vicinity of 
the hole satisfies the boundary conditions 

lim \dt + iCm/D^ - drA i'm = 0, (10) 

lim [dt+dr,]^rn = 0. (11) 
r.^ — f +00 



B. Geodesies 



Hence, by investigating the null geodesies of the line ele- 
ment ([5]) we may understand high-frequency wave prop- 
agation. Let us consider a geodesic with tangent vector 
k'^ = (i,f, (/)), where the overdot denotes differentiation 
with respect to an affine parameter. Geodesies have two 
constants of motion, i.e., energy and angular momentum 



L = -Ci + r^(j), 

respectively. For a null geodesic {k^kf^ — 0) we may write 
the 'energy equation' as 



r^ + Vg,o{r)^E^, (13) 



where 



^gco(0 = ( 1 - — ^ + — ^. (14) 



Let us consider a null geodesic impinging from infinity. 
Its 'impact parameter' b may be defined as the perpendic- 
ular distance (measured at infinity) between the geodesic 
and a parallel line that passes through the origin, 

b=L/E + C. 

Here b is defined as a displacement that may take ei- 
ther sign, allowing us to distinguish between co-rotating 
{+) and counter-rotating (— ) geodesies. Now, if |&| is 
large, the geodesic will be scattered; if |6| is small, the 
geodesic will be absorbed (i.e. it will pass through the 
horizon). Between these regimes is a 'critical' geodesic, 
which is neither scattered nor absorbed but instead ends 
in perpetual orbit at r = r^. By solving simultaneous 
conditions V^co(?'c) = E'^ and dVgco/dr{rc) — we find a 
pair of 'critical' geodesies given by 



+ C2 \bf -C\ 



1/2 



(15) 
(16) 



Note that b^ is defined to be negative. For C > 0, 
geodesies with b > C {L > 0) co-rotate with the system, 
whereas geodesies with b < C (L < 0) counter-rotate. 

It is natural to define I = b~C, which is positive for co- 
rotating geodesies and negative for the counter-rotating 
ones, so that: 



l^=bf~C = ~2C ± 2v/l>2 + C'2. 



(17) 



Trajectories of null geodesies impinging from spatial in- 
finity upon a draining bathtub are illustrated in Ref . |25j . 



According to the eikonal approximation, very short- 
wavelength perturbations propagate along null geodesies 
of the effective spacetime ([5]), i.e. 



$ ~ exp i-ik^xf") , fc^fc^ = 0, fc^fc^. 



(12) 



III. QUASINORMAL MODE RESONANCES 

With these preliminaries established, let us now turn 
our attention to quasinormal (QN) modes. We start in 
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Sec. Ill A by demonstrating QN resonances in a perturba- 



tion encroaching upon a black hole, by evolving Gaussian 
initial data in the time domain. In Sec. |IIIB| we recap 
the theory of QN modes as poles of the Green function 
in the frequency domain, and in Sec. Ill C we describe 
the symmetries of the spectrum and provide an exact so- 
lution for the non-resonant m = mode. Next, we apply 
two frequency-domain methods to determine the spec- 
trum: the numerical method of [51], and the 'geodesic 
expansion' method of [29ll3Q]. We validate our results, 
compare with the Kerr spectrum, and give a geometric 
interpretation. Henceforth we set = D = 1 for conve- 
nience, unless otherwise stated. 



^ 0.0001 - 



Simulation of quasinormai mode ringing in a draining bathtub 

mode m - 3 




A. Quasinormai modes in the time domain 

To study the phenomenon of QN mode ringing, let 
us consider the evolution of a small perturbation to the 
system. For concreteness, we take an initial condition of 
the form 



ipmit = 0,r*) 
dt-tpmit = 0,r) 



exp[-(r, -r,o)V(2<T2)], (18) 
0, (19) 



where r^,g and cr^ are arbitrary constants. Now, we sim- 
ulate the evolution of this initial perturbation by apply- 
ing finite-difference (FD) methods to the wave equation 
^. There are many possible FD methods; we chose the 
'Method of Lines' (described e.g. in Ref. and Sec. 4.2 
of Ref. [32]), using second-order differencing on spatial 
slices and the fourth-order Runge-Kutta method to ad- 
vance in time. This method was chosen primarily for its 
good stability properties: it is numerically stable pro- 
vided that the time step r is small enough (typically we 
took T = h/2 for our evolutions, where h is the spatial 
grid spacing) . We made the spatial domain large enough 
in so that the 'physical' perturbation did not encounter 
the spatial boundaries during the simulation run. 

Typical results are given in Fig. [l] which shows the 
reponse il>m{t,T) at fixed radius r — lOr/j as a function 
of time, for an initial condition of the form (19). The 



logarithmic scale on the vertical axis shows that, at in- 
termediate times, the response is apparently dominated 
by exponentially-damped ringing. This response is typi- 
cal of systems with an unstable geodesic orbit (or equiv- 
alently a peak in the potential barrier). Ringing may 
be understood in terms of QN modes with characteristic 
complex frequencies which depends on parameters 
of the system (C and D) rather than details of the initial 
perturbation. 

The time-domain signal in Fig. [l] clearly shows the im- 
print of the least-damped QN mode. In later sections, we 
find that (for m ^ 0) the least-damped mode is the ret- 
rograde (rather than prograde) mode. Fig. [T] shows that 
the damping rate and ringing frequency of this mode de- 
creases as the circulation rate C increases; in Sec. jlll E| 
we relate this behavior to the properties of geodesies. 



FIG. 1: Illustrating the 'quasinormai ringing' response as a 
function of time t of a draining bathtub to a Gaussian initial 
perturbation. Here the value of il}{t, r) at a fixed radius r = 
lOr^ is shown, for various values of the 'circulation' C of the 
bathtub, C = 0...2. As the circulation rate increases, the 
ringing frequency and the damping rate decrease. 



The ringing frequency increases with mode number m. 
Figure [2] shows that, at very late times, the signal is 
dominated by a power law decay. 



V'(t,r) 



where 




(20) 



In other words, the decay is one power of t more rapid in 
the special case of time-symmetric initial data [such as 
Eq. (19)]. For further discussion of power-law decay in 
the DBT context, see Sec. HE in Ref. [20j . 



B. Quasinormai modes in the frequency domain 

It is well-known that the origin of QN ringing and 
power law decay may be understood by considering the 
decomposition of the Green's function in the frequency 
domain (see for example [33l [34]). Let us briefly re- 
cap the argument here. A formal solution of Eq. ^\ 
with arbitrary initial condition '0^(a;) = '0m(^ = 0,xj, 
'il)%Xx) = dt'ipmit, x)\t=o may be written as 



G(x,a;',tX(a;')+ 



dtG{x,x',t) tPl^ {x')] dx', 



(21) 



where x = r^, x' — r'^^ and G is the (retarded) Green's 
function defined by 



/ d iCm 



VW(r) 



G^6{t)S{x-x'), 
(22) 



5 



Power law decay for modes of a draining bathtub 
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FIG. 2: Log-log plot illustrating 'power law decay' in modes 
m — ... 3 for a perturbation of a draining bathtub. At 
late times, the perturbation decays as t/i cx {t/rh)~^, where 
rj = 2\m\ + 1 for generic (in this case, Gaussian) initial data, 
and ?7 = 2|m| + 2 for time-symmetric initial data, dttp{t = 
0, r) = 0, which is shown here. The decay rate is independent 
of the circulation C. 



and G = for t < 0. The Green's function may be 
written in terms of an inverse Fourier transform 

G{x,x',t) = -- Gix,x',u;)e-''^'duj, (23) 

./-oo+iC 



with C positive and real. The frequency-domain function 
G is constructed in the following way: 



G{x, x ,uj) 



1 <^(^XPc.(^0, ^<^', 

'2ii^AZ^\i^^^{x')u-^^{x), x>x'. 

(24) 

Here u™^{r) ('in') and u^^^ ('up') are solutions of the 
homogeneous equation 



Cm 



u{r) = 



(25) 



that are, respectively, ingoing at the horizon and outgo- 
ing at infinity, satisfying the following asymptotic bound- 
ary conditions 



and 



r>ont ^-\-iCjr^ _i_ Din „ — icjr* 



(26) 

r* — >■ — oo, 
r* — >• -t-oo, 

(27) 

where w = a; — mC jD"^ . By considering the Wronskian, 
one may establish certain relationships between these 
constants, for example, wA'^^ = ujB'^l^. 

T he f requency-domain Green's function G given in 
Eq. (24 1 has poles at frequencies for which = 0. 



Such poles do not lie on the real frequency axis, but 
rather in the lower half of the complex frequency plane. 
For t > 0, the contour of integration in Eq. (23) may 
be closed in the lower half-plane, as described e.g. in 
Refs. [331131], enclosing the poles. By Cauchy's theorem 
there arises a sum over residues of these poles, which 
is known as a 'QN mode sum'. The QN mode sum 
manifests itself as the damped ringing response seen e.g. 
in Fig. [l] Furthermore, there exists a branch point in G 
at zero frequency, which necessitates a branch cut along 
the negative imaginary axis. The integral of G along ei- 
ther side of the branch cut is associated with the power- 
law decay at late times, as observed in Fig. [2j With this 
in mind, let us now consider the QN spectrum, i.e. the 
set of frequencies defined by AJ"^ — 0. 



C. The quasinormal mode spectrum 

1. The m = mode 

The TO = mode is isotropic and independent of C, 
and furthermore, Eq. ( 25 ) has a simple closed-form solu- 



tion in this case. The 'in' mode is 

<M = rl/2e-'^(^^)»/^^(l - */3)/_,^(*c.r/i/2), (28) 

where /3 = loD and Iv{z) is a modified Bessel function of 
the first kind [35j . In this case, the boundary condition 
constants are 



= (27rz/3)-^/2e-^'^(*/3)'^^r(l-*/3), 



*6 ■ 



(29) 
(30) 



The 'up' mode is 

up^_^^ _ (2z7r^r) 

2 sinh(7r/3) 



"oS('') = 



e^^/+.^(zc.r/i/2) 



(31) 



From Eq. (30) we infer that A^§^^ has no zeros in the 
complex-cj plane. This leads us to conclude that there 
are no quasinormal modes for m = 0. The solutions for 
higher modes to ^ also have closed form expressions, 
but in terms of the lesser-known confluent Heun functions 
(see e.g. [36]). 



2. Symmetries of the spectrum 

Through our time-domain simulations we found that 
QN ringing is a key feature in the response of higher 
modes (to ^ 0). Although only one QN mode fre- 
quency is apparent in the time-domain data of Fig. [T] 
for each mode m ^ there is in fact an infinite number 
of damped overtones present in the QN spectrum, labeled 
by n = 0, 1, . . . , oo. For a given m, n, there is a pair of 
modes: one co-rotating (-I-) and one counter-rotating (— ) 
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with the circulating flow, with frequencies a;^„(C). The 
spectrum has the foUowing symmetries: 

a.±,„(C) = ~^^*^,AC) (32) 

and 

^±,„(C)=c.f,„,„(-C), (33) 

where * denotes complex-conjugation. Note that (for 
C > 0) the co-rotating mode oscillates and decays 
faster than the counter-rotating mode, i.e. |Re(w^„)| > 
|R-e(w-„)| and |Im(a;+„)| > |Im(a;-„)| (with equality if 
C = 0), and the frequencies scale linearly with 1/rh- 



the imaginary part of the least-damped (n = 0) counter- 
rotating mode asymptotes to zero, whereas the imaginary 
part of the co-rotating mode increases without bound. 

Somewhat different behavior is observed in the spec- 
trum of equatorial modes (m ~ I) of the scalar field of 
the Kerr black hole (BH), illustrated in the right plot of 
Fig. [3] In the Kerr case, the relevant rotation parameter 
is a = J/M (rather than C), where J and M are the 
angular momentum and mass of the BH, respectively. 
As in the DBT case, the co-rotating (-1-) modes oscil- 
late more rapidly than the counter-rotating modes (— ), 
|Re(u;+„„)| > \R.e{LOi^J\. However, unlike the DBT, the 
imaginary part of both modes (±) is found to decrease 
in magnitude as the rotation rate increases. 



D. The continued-fraction method 

The continued fraction method is a fast and accurate 
numerical method for determining QN frequencies. It 
was first applied to determine the QN spectrum of black 
holes in Ref. [37] , and was adapted to the draining bath- 
tub case in Ref. [21]. Here we briefly recap the method 
giving some results in Sec. |III G[ 

The starting point is an ansatz for the QN wavefunc- 
tion, 



r- 1 
r -M 



-iu/2 oo 



k=0 



(34) 



It was shown in Ref. [H] that the coefficients Ok satisfy 
a four-term recurrence relation, and, by making use of 
Gaussian elimination, it may be reduced to a three-term 
relation, 



a/cOfc+i -I- /^fcflfc -I- 7fcafc- 



0. 



(35) 



Here ak, Pk and 7fc, are complex coefficients that depend 
upon the frequency lo and also upon m, C and D. For 
Eq. (|34| to represent a valid QN mode solution, the sum 
Ok must converge to a finite value. This condition is 
equivalent to the following continued-fraction condition: 



/3o 



a-oli ai72 ^273 
A- /32- ft- 



= 0. 



(36) 



Finding the roots of Eq. ( 36 ) in the complex-frequency 



domain is a straightforward task for a numerical min- 
imization algorithm. High numerical accuracy for the 
QN frequencies may be obtained. 

In the left plot of Fig. |3] we show the dependence of 
the 'fundamental' (n = 0) mode on m and on the rota- 
tion rate C of the acoustic hole. The spectrum exhibits 
the symmetries (32) and (33 1. In the non-rotating case 
(C = 0) the two modes and a;~„ are symmetric 

under reflection in the imaginary axis. The introduc- 
tion of rotation (C ^ 0) breaks this symmetry. The co- 
rotating (-I-) mode increases the magnitude of its real and 
imaginary parts, whereas the counter-rotating mode (— ) 
moves in the opposite way. In the limit of very large C, 



E. Geometric interpretation 

This behaviour may be understood through an approx- 
imate formula for the QN frequencies [38l [39] , 



n^m-iK^{n+ 1/2). 



(37) 



Here, is the orbital frequency of the prograde (-I-) or 
retrograde (— ) null orbit, and is the corresponding 
Lyapunov exponent |40| I41j. In the case of the DBT, we 
have 



I /If, 



A' 



1/r, 



± 

c ' 



(38) 



where and rf are given in Eqs. (17) and (16), respec- 
tively. Since r~ and /~ increase in magnitude with C, 
the magnitude of the real and imaginary parts of a;~„ 
decreases with C (and the opposite is true for r+ , and 
i-^mn)- By comparison, for the equatorial modes of the 
Kerr BH one has instead 



where 



= 1/5; 



2M 



C ' 



l~2a/bi 



A± = 



cos I - COS ^(=Fa/A/) 



bf = ±3\/Mrt - a, 
for the equatorial orbit [39 . Expanding A* gives 



A± = (27) 



-1/2 



^ 2a2 10^3 3 , 4, 



(39) 

(40) 
(41) 

(42) 



Clearly, for both co- and counter-rotating equatorial or- 
bits on Kerr, A* decreases with \a\. This 'explains' the 
observation in the right plot of Fig. [3] that the damping 
decreases with \a\. 

In Ref. [30| . a more accurate approximation for the 
QN modes of the Kerr BH was found using a 'geodesic 
expansion method', which builds upon the understanding 
of the properties of the null orbits. Let us now apply this 
method in the DBT case. 
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Quasinormal mode spectrum of the draining bathitub 



Quasinorma! mode spectrum of equatorial {[m|=l) modes of ttie scalar field or Kerr spacefime 



C = 
C = 0.2 • 



C = 
C = 0.2 



-1 1 

Re(cor^) 




-2 -1.5 



FIG. 3: Left: QN mode spectrum of the draining bathtub. Right: QN mode spectrum of the equatorial modes of the scalar 
field on Kerr spacetime. The left plot shows the fundamental (n = 0) quasinormal mode frequencies io^„ of modes m = ±1 
and ±2, for a range of circulation rates, C = 0, 0.2, 0.4, . . . 2.0. For each mode m 7^ there are two QN frequencies uj^. 
In the non-circulating case (C = 0), the spectrum has the symmetry = — w^*- Circulation creates a difference between 
co-rotating (a;+) and counter-rotating (oi^) modes: Counter-rotating modes oscillate and decay more slowly t han co-rotating 
modes (|Re(aj+)| > |Re(a;~)| and |Im(a;^)| > |Im(ctj^)|). The spectrum obeys the symmetries of Eqs. (321 and (33 1. The right 
plot shows the fundamental (n = 0) quasinormal mode frequencies to^^rh of equatorial (|m| = I) modes m = 0, ±1 and ±2 of 
the scalar field on the Kerr spacetime, for a range of the black hole rotation rates a/M = 0,0.2,0.4,0.6,0.8,0.99. Co-rotating 
modes oscillate more rapidly than counter-rotating modes, |Re(cjj^„)| > |Re(cj;~^„)|. In general, both co- and counter-rotating 
modes become less damped as the rotation rate increases. 



F. The geodesic expansion method 

First let us rewrite the function u (r) appearing in 
Eq. ( 25 ) using the following ansatz [53 [10] ■ 



Urn (r) = X (r) exp 



where 



a = w — 



Cr 



1 - 



CL 



1 - 



(43) 



(44) 



Note that here 



= and L 



- ;± 



, i.e. representing 



either the co- or counter-rotating cases, as defined in Eqs. 
(16) and (17), respectively; we drop the ± superscript 



here for clarity. Substituting ansatz (43) into the radial 



equation ( 25 ) leads to an equation for the function x {'"') j 
namely 



fx" + (/' + 2ia) x' + 




(45) 



where ' denotes differentiation with respect to r. Now, 
to seek the QN frequencies LOmn and radial wave func- 



tions X. (^) we introduce an expansion in terms of inverse 
powers of m, namely 



X = 



9=-l 

n 00 

00 

x[]exp(m-«5(")(r)), 
9=0 



where 1^ = 1 — r*/r. Here vjq'"' and a|"^ are dimension- 
less coefficients, and Sq'^\r) are smooth radial functions. 
To determine these unknowns, we impose a condition of 
regularity upon the solution at r ^ Vc- 

Let us illustrate the approach by computing the QN 
frequencies and radial wave functions for the fundamen- 
tal mode (n = 0) for the co-rotating case. Substituting 
the expansions in Eqs. (46) and (47) into the radial equa- 



(46) 



(47) 



tion ( 45 ) , and then rewriting it in terms of the powers of 



m associated with coefficients of the expansion, leads to 
the following system of equations: 
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i^iizzz^zL _ 1 = 

2rl{r^-[rl-ll/2]) 2» / 

IcT'' lc\ 



1, 



1 - 



= 0, 



(48) 
(49) 



r 



etc., where the superscript n has been suppressed for 
simphcity. Equation (48) is identically satisfied through 
the choice w^i ~ 1. 
S'o(r) is regular at r 



To solve Eq. 



Tc, which leads to vjq 
and to the first-order differential equation 



(49) we impose that 



(0) 



(r-rc)(r + re)^ 
ar 



r 
' c 

r 



'■c 

2rr 



1 - 



(50) 



The radial function may be obtained by integrating 
this equation. In subsequent equations at orders ©(m") 



and higher, the second derivative Sq appears, which may 
be obtained by differentiating (50). The higher-order 



equations are solved in a similar way: taking the equa- 
tion at 0{m^~'') we first impose the continuity condition 
at r = Tc to obtain 137^, and next solve to obtain S"^. 

Using a symbolic algebra package (e.g. Maple or Math- 
ematica), this procedure may be automated, and the ex- 
pansion may be taken to higher orders. We have com- 
puted the expansion of the QN frequency up to the order 
0(rn~^). Below, we quote the expansion of the frequency 
of a general overtone n (= N —1/2) up to order 0(m~''): 



In OJ' 



ilr 



c ^mn 



N + mr^ 



Tc Ll28r4 
^ ll{2ill-SQrl)^N^ 



4 V c c) 32 r4 



1024 



4- m 



1 



1048576 



4096 



ll{bli~imlrl + ?,Mrl 



(64640/2r^ - 21040r^/;[ + 2\2bll - 57344r^) 



^^It {Mmlrt ~ QlQrllt + 75/6 - 4096r6) + ^895^4 _ 592^2^2 ^ 91^4^ ^4 



(51) 



Note that, for a given m, this expression yields two 
QN frequencies: obtained using co-rotating geodesic 
parameters {lt,^t) ^'^d uj^n obtained using counter- 
rotating geodesic parameters {l~ , r~). The frequency ex- 
pansion exhibits the symmetries ( |32[ ) and ( 33 ) , as may be 
confirmed with the aid of the relations r^JC) — rj(— C) 
and lf{C) = -l^{-C). 



G. Validation 

In Table |T] we compare the QN frequencies found via 
Eq. (51) with numerically-accurate values obtained via 
the continued-fraction method of Ref. [HI], and re- 
obtained here in Sec. |III D[ At large m, we find excel- 



lent agreement. At small m, the approximation is not so 
good. For instance, for m = 1 the most accurate estimate 
was found by truncating the series at 0(m~^), suggest- 
ing that Eq. (51) is in fact an asymptotic series. In 
Fig. [4] we plot the difference between continued-fraction 
and geodesic-expansion results, to confirm that Eq. (51) 



is indeed valid to the stated order in the large-|m| regime. 



IV. REGGE POLE RESONANCES 

Regge poles (RPs) are closely related to QN modes 
[42] . Both types of resonance are associated with the 
zeros of A™^^ . Whereas QN modes occur at real values of 
m (and complex oj), the Regge poles (RPs) occur for real 
values of w (and complex m). That is, m^n is a Regge 
pole angular momentum if 



AZui.^-.'n^ojn) = 0. 



(52) 



Methods 



Given the close relationship between QN modes and 
RPs, it is no surprise to find that methods used for QN 
modes can be easily adapted to locate RPs ^ 



For instance, the continued-fraction method (Sec. HID) 
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TABLE I: QN frequencies of the fundamental mode (n = 0) for m = 1 . . . 5. The second column gives the frequencies determined 
via the continued-fraction method of Sec. |IIID| The third column gives the frequencies estimated from the expansion method, 
Eq. (51 1. The numeral in parentheses indicates the absolute error in the last displayed digit, where the error in the expansion 
method was estimated from the magnitude of the final terms in the series, Eq. (51 1. Note that the value marked with an 
asterisk (*) was obtained by truncating the series at 0{m~^). 



D = 1, C — 0, symmetric modes (±) 



m 

1 

2 

3 

4 

5 



Continued- Fraction 



Expansion Method 



±0.4068326196672(2) - 
±0.95272808772474(6) 
±1.46854069662523(5) 
±1.9764527143560(1) - 
±2.4811874975616(1) - 



0.341236118125(2)i 

- 0.3507394957317(2)i 

- 0.3524255329360(2)i 
0.3529594212624(2)j 
0.35318833783278(6)i 



±0.40(1)* - 0.345(8)i* 
±0.9524(9) - 0.3511(2)i 
±1.46854(6) - 0.352444(8)i 
±1.976453(8) - 0.352961(l)i 
±2.481188(2) - 0.3531886(l)i 



D — 1, C — 1, counter-rotating modes (— ) 



m Continued-Fraction 



Expansion Method 



1 -0.1490148506555(1) 

2 -0.3875750711220(2) 

3 -0.6040167016214(2) 

4 -0.8155812868761(1) 

5 -1.0253078467593(1) 



0.1945123420127(2)i 

0.1930307208974(l)i 

0.19220554998274(l)j 

0.1918516723182(3)i 

0.1916753960108(2)1 



-0.14(1) - 0.200(8)1 
-0.38753(6) - 0.19296(9)i 
-0.604014(4) - 0.192204(3)1 
-0.8155811(5) - 0.1918515(3)i 
-1.0253078(1) - 0.19167538(6)i 



D = 1, C = 1, co-rotating modes (±) 



m 

1 

2 

3 

4 

5 



Continued- Fraction 



Expansion Method 



±1.13097639081141(2) 
±2.3742286767072(2) - 
±3.5943282234311(3) - 
±4.8080846946733(2) - 
±6.0192217565670(3) - 



- 0.4485344193663(3)1 

0.45797306346650(6)1 

0.4600949990202(3)i 

0.46088327915183(4)1 

0.4612576561145(2)i 



±1.134(9) - 0.445(7)j 
±2.37421(3) - 0.45791(2)1 
±3.594327(2) - 0.460093(2)i 
±4.8080845(3) - 0.4608831(2)i 
±6.01922173(6) - 0.46125763(3)j 



Validation of tine frequency expansion 






|Re(a) 

|im(a) 

IRejA) 

llmja) 


.C-0 
,C.O X 
,C = 1 * 
,C.1 □ 



































may be used without modification, if m is allowed to 
become complex. Furthermore, the geodesic expansion 



method (Sec. IIIF) can be easily modified to find RPs 
[29] . Either one may repeat the arguments of Sec. IIIF 



or one may take the following simple steps: (i) assume 
that the Regge poles have the expansion of the form 



C-iw±Co±Ci^ + 



(53) 



FIG. 4: Validation of the frequency expansion, Eq. (51 1. The 
log-scaled plot shows the real and imaginary parts of A as a 
function of mode number m, for the cases C = and C = 1 



(and D = 1). Here A : 



; LUCF - 



where t^cF is the QN fre- 



quency obtained via the continued fraction method and LOc^p 
is the frequency estimate given by Eq. (51 1. The plot shows 
that the imaginary (real) part of A scales as m"* (m~^), as 
predicted by Eq. (51 1. 



(ii) substitute the expansion for QN frequencies into the 
above equation, and (iii) solve order-by-order in m to find 
the expansion coefficients (k- We find 
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lG24rJ 



32r4 



128r4 



iN 



Ic 

4096rZ 



Ic (-73728r6 + 74880r4/2 _ 22640r2/4 + 2125/^ 



Ic (36864r6 - 7808r^/2 _ IGlGr^/^ + 75^6 



131072r. 



10 



1048576rio 
ll [iMrl + bmrlll + 455/^) 



65536r. 
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+ 0(^-4)^ (54) 



1/2. Again, as for Eq. (51), we note that 
— ± 



obtained 



where N = n 

this expression yields two separate values 
by using either co-rotating or counter-rotating geodesic 
parameters (/^,r^). We note that it is straightforward 
to take the expansion to much higher orders if desired. 
Finally, it is relatively straightforward to find Regge poles 
by direct integration methods, because (unlike quasinor- 
mal modes), RP modes are not divergent in the limits 
r* — )■ ±c». 



B. Spectrum 

The Regge pole spectrum has the following symmetries 

mlJC) = -m±:^„(C) (55) 

and 



3 pole spectrum of draining bathtub 



(56) 



The Regge pole spectrum in the complex-m plane is 
illustrated in Fig. [5] For w > 0, the co-rotating (-I-) 
modes lie in the first quadrant, and the counter- rotating 
(— ) modes lie in the fourth quadrant. In the non-rotating 
case, the spectrum is symmetric, m+„ = —m~„. For 
C > 0, the real and imaginary parts of r7T,+„ decrease 
with C, whereas the real and imaginary parts of — ?tiJ„ 
increase with C, as expected from Eq. (54). 



V. APPLICATIONS OF COMPLEX ANGULAR 
MOMENTUM (CAM) METHOD 

In this section we show that Regge poles play an 
important role in both absorption and scattering pro- 
cesses; their role is revealed through application of the 
so-called Complex Angular Momentum (CAM) method. 
The CAM method makes a link between certain oscilla- 
tions in cross sections (see below) and the poles of the 
scattering matrix S{m) in the complex-m plane, i.e. the 
Regge poles. The CAM method was successfully applied 
to scattering by a Schwarzschild black hole in [23], and 
was recently extended to treat absorption [27, 28, 46 . We 
believe this represents its first application to a rotating 
spacetime. 
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FIG. 5: Regge pole spectrum of the draining bathtub. For 
three frequencies \iu\ — 0.25, 0.5, and 1.0, we plot points for 
rotation rates C/r^ = 0, 0.1, 0.2, . . . , 1, showing the funda- 
mental (n = 0) Regge pole values m^^r^ in the complex-m 
plane. Co-rotating and coun ter- rotati ng modes are shown; 
they exhibit the symmetries (55) and (56). For C > 0, the 



co-rotating (counter-rotating) modes lie in the upper (lower) 
half plane. 



A. Absorption Cross Section 

The absorption cross section (Tabs of the draining bath- 
tub was studied in Ref. |25j . It may be obtained from a 
sum over partial wave contributions, via 



Cabs — 



00 

- E 



m— — oo 



(57) 



At high frequencies, the absorption cross section is 
seen to approach the 'geometric' capture cross section, 
CT^™) = \h+\ + \h-\ = 4VL>2 + C2. In the non-rotating 
case, C = 0, (Tabs displays regular damped oscillations 
with increasing frequency. With increasing C, these os- 
cillations become less regular (cf. Fig. [t]). 

Recent papers [171 [551 US] on the universality of 
high-frequency absorption cross sections for spherically- 
symmetric black holes use CAM methods [22l [23l l42l |47] 
to show how oscillations in the absorption cross section 
(with frequency uj) are directly related to the Regge pole 
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spectrum. 

To apply CAM methods, we must employ a 
analytic continuation into the complex-m plane 
therefore define 


suitable 
Let us 




lm(m) 


♦ 

.^n=l 


( + ) 


OO 

o-abs = - V r(m), 


(58) 


C+ 




^n=0 




> 






in— — OQ 










where 











C- 


ylOUt / /lout \* 

Tim) — 1 mujy m-uj) 

4'" ('4'" \*' 


(59) 


(-) * 









-Re(m) 



A series such as the one above can be expressed as a 
contour integral using the Watson transformation, i.e. 



fabs 



2a; Jq sin(7rTO) 



r(m)(im. 



(60) 



where now m takes complex values. Here the contour of 
integration C encloses in a clockwise sense all (and only) 
the poles of the integrand that lie on integer values on 
the real axis at m = — oo, . . . , +cxd, i.e. 

/• />oo±ic 

, where / — / ' (61) 

C JC+ JC- Jc± J-oo±ic 

for some small positive real value c. See Fig. [6] for a 
graphical representation of the contour. Next, we use 
the identity e"™ — e""^'^ + 2i sin(7rTO) to obtain 



(int) 
0-abs = (Tabs 



T(to) 



2uj Jq sin(7rm) 



dm 



^r(m) 



2uj Jq sin(7rm) 



dm. 



(62) 



where 



(int) 



^abs - - / r(m)dm. (63) 

UJ 



The contour C-^- may be closed in the upper half-plane, 
enclosing poles of r(m) at m+„ and m~*. Likewise, the 
contour C_ may be closed in the lower half-plane, enclos- 
ing the poles of r(m) at mj„ and m^*. This leads to a 
sum of residues, that is. 



Cabs — CTabs 



2-ReEE ' 



± n-. 



^ sin [-Kmtn) 



Here 7^„ denotes the residue of r(TO), 
7±,- Res^ [r(m)] 



(64) 



(65) 



and we have made use of the symmetry relation 

Res^ [r(m)] = (7±„)*. (66) 



FIG. 6: Contour integration in the complex-m plane. The 
circles show the positions of the poles of 1/ sin(m7r), at integer 
values of m. The stars illustrate the positions of the poles (for 
> 0, C > 0) of S{m), i.e. the Regge poles, for overtones 
n — . . . <x. The contour C+ is closed in the upper half-plane, 
and the contour C_ is closed in the lower half-plane. 



Equation ( 64 ) is an exact expression which may be com- 
puted numerically. The sum over overtones converges 
rapidly, due to the damping effect of the imaginary part 
of the Regge pole. The residues may be computed with 
a numerical scheme based on integration of the radial 
equation (25), using 



' ( A° 



(67) 



where A^^^ « a^„(m — m^^) + . . . in the vicinity of a 
Regge pole. The integral over frequency can also be per- 
formed in a straightforward manner. Numerical results 
are shown in Fig. [7] We find a very good agreement with 
the results of the partial wave method of Ref. ^25] . 

To better understand the geometric meaning of 
Eq. (64) in the high-frequency limit, we may proceed by 



making some simplistic approximations: 

r(m) Q{l'^u! — m)Q{m — 1^^), 
where 0(-) is the Heaviside step function, and 



and furthermore 



27r 



(68) 



(69) 



(70) 



where r]"^ — Neglecting the higher overtones 

leads to the simple approximation 



abs 



l^\ l-h47r?7±e-''''*sinc(27r|/±|a;) 



(71) 



wherejjabs = f^ths + ''^abs ^^'-^ sinc(a;) = sin(x)/a;. Equa- 



tion ( 71 ) makes it clear that the cross section approaches 
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the geometric capture cross section in the high-frequency 
hmit. Furthermore, the cross section exhibits damped 
oscillations with increasing frequency, which are due to 
the Regge pole contribution, which are controlled by the 
geometric quantities {l^ and r^). 

Note that a somewhat more accurate approximation 
may be obtained by using an improved approximation 
(cf. Eq. (19) in Ref. [27') for the transmission factor. 



r(m) 



1 + exp 1^— 27r 
1 + exp ( — 27r 



2mr/+ 

[Me')'-'"' 



m > 0, 
TO < 0. 



(72) 

In Fig. [7]we compare the results of the partial- wave cal- 
culation of Ref. [25] with a numerical calculation based 
on the exact Regge pole expression, Eq. ( [64| , and the 
'geometric' high-frequency approximation, Eq. (71 1. It 
is clear from the lower right plot of Fig. [7] that the co- 
rotating (counter-rotating) orbits are responsible for os- 
cillations in the profile with a higher (lower) frequency 
and smaller (larger) amplitude. The superposition of 
CO- and counter-rotating effects creates a rather irregular 
profile, in contrast with the more regular profile exhibited 
in the non- rotating case (left plots). 



B. Scattering cross sections and 'orbiting' 

Various studies of monochromatic planar-wave scat- 
tering by black holes [15H5T] have shown that there may 
arise regular oscillations in the scattering cross section, 
whose angular width is inversely proportional to incident 
wavelength. This phenomenon is called 'orbiting' [52^i53]. 
or alternatively, 'spiral scattering' [35]. In essence, such 
oscillations arise due to interference between pairs of rays 
which pass in opposite senses around the black hole. A 
study of scattering by the DBT [SH [55] has recently 
demonstrated that orbiting oscillations also arise in this 
case. 

(Note that, in three-dimensional scattering, orbiting 
is supplemented by another effect: the 'glory'. A glory 
is a bright spot (or ring) in the forward- or backward- 
scattering directions, whose width (amplitude) decreases 
(increases) linearly with frequency. Semi-classically, a 
glory is created by a one-parameter family of geodesies 
scattering into the same solid angle. This is not possible 
in 2D; hence the glory effect is absent in the DBT case.) 

In this section, we apply the CAM method to obtain a 
deeper understanding of orbiting oscillations in the scat- 
tering cross section. Let us begin with the partial-wave 
expression for the scattering amplitude, 



m— — oo 



[S{r 



1]' 



(73) 



where k — {2iTTUj) and the 'scattering matrix' S'(m) 
(which is scalar-valued in this case) is defined by 



Sim) ^ 



(74) 



We note that, with this definition, the scattering matrix 
has the following symmetry: 



S{^m, C) = e-2'""5'(TO, -C). 



(75) 



The differential cross section is simply da/dcj) = 

We will adopt the convention that the scattering angle is 

in the range < < 27r (rather than — tt < < tt). 

As in Sec. |V A[ the key step is to apply a Watson trans- 
formation 52 to convert the partial wave series ( 73 ) into 
a contour integral: 



^ in f 5(TO)e™'^ 

2 Jq sm(7rTO) 



dm, 



(76) 



where ^(to) = e-"™5(TO) 



Let us proceed by recasting ( 76 1 into the form 

f -t-oo 



IK, 



S{m)e''^^*-'^'''^dm 
S{m)e'"^'^ 



sin(7rTO) 



-dm — 



IK 

~2 



S{m)e 



zm(0— 27r) 



sin(7rm) 



-dm. 



In the high-frequency limit, we may neglect the first in- 
tegral, because it has no stationary phase points in the 
range < < 27r. Furthermore, we may evaluate the 
second and third terms by closing the contours C+ and 
C_ in the upper and lower half-planes, respectively. Ap- 
plying Cauchy's theorem leads us to 



(77) 



su\[-KmZn ) sm(7rTOaj„ ) 

where the residues are defined as follows. 

Res S{m). (78) 



The 'fundamental' mode {n = 0) gives the dominant con- 
tribution, and the sum over overtones appears to be ex- 
ponentially convergent. 



1. High-frequency approximation 

The geodesic expansion method led us to expressions 
for the Regge poles and corresponding wave functions, 
which are valid in the high-frequency regime. Now, by 
combining this method with standard WKB techniques, 
we may also obtain a low-order approximation for the 
residues (78 1. As the calculation is rather lengthy, we 



quote the key results here, and give a fuller exposition in 
Appendix [A] 

In the high-frequency limit, the RP wave function is 
approximately 





-rc\ 


r 


( r \ 











exp [iujr — iib(r^ — r)) , 
(79) 
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Absorption cross section of non-rotating bathtub 



Absorption cross section of rotating batfitub, C = 1 



Partiai wave metiiod 
Compiex anguiar momentum method 
Geometric Approximation 
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FIG. 7: Absorption cross section as a function of the frequency, for C = (left plots) and C jru ~ 1 (right plots). In the upper 
plots, the solid line shows the absorption cross section crabs computed from the partial wave sum [25], the dashed line shows 
o"abs computed numerically from the complex angular momentum expression [Eq. ( 64 1] . Note that these two lines lie very close 
together. The dotted line shows the simple geometric approximation [Eq. \71} ] , which is valid in the high-frequency regime. The 
lower plots show the separate contributions from Regge pole sums and the co-rotating and counter-rotating integrals [Eq. ( 62 1] . 



where uj 



Here we have chosen the nor- 



malisation such that Uuynif) ~ e*' 
The residue is approximately 



in the large-r limit. 



±8 (-iwrc) J 



where 



A = 



ijj — In 
2 \rr 



(27r)i/2r 



rh 

Th 



(80) 



(81) 



Orbiting oscillations in the cross section |/cj|^ arise 
from interference between co-rotating (+) and counter- 
rotating (— ) contributions to (77 1. The angular oscil- 
lation frequency is given by Re (m^„ — mj„) « i^(^^ + 
l^^l) — Aujre at lowest order (where = y/C^ + D'^/c). 
Therefore, the angular width of the orbiting oscillations 
at high-frequency is simply 



2ujr„ 



(82) 



The imaginary part of m+„ (or m^„) determines the rate 
of attenuation with scattering angle (or 2tt — (p). 



2. Numerical results 



Figure [8] shows a comparison between the high- 
frequency approximation for the residues ( 80 ) , and ac- 



curate values computed using a numerical method based 
on direct integration. We find a good agreement, sug- 
gesting that the analytic approximation Eq. ( 80 ) can be 
used to capture the key features of orbiting. 

Figure [9] compares the CAM approximation for orbit- 
ing ( |77[ ), with accurate numerical results obtained by 
summing the partial wave series [SS]. Two CAM ap- 
proximations are shown; for the first approximation we 
computed the RPs and residues numerically (for low over- 
tones n = ... 4) , and for the second approximation we 



used the asymptotic analytic results, (54) and (80) 



The agreement is found to be good, particularly at 
large scattering angles. In other words, the Regge pole 
approximation, given in ( 77 ) , accounts rather well for the 



orbiting oscillations. In particular, the asymptotic results 
(given in closed form in terms of geodesic parameters 
and l^) provide an excellent description in the semi- 
classical (high-frequency) regime. 
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Residues of Regge poles of draining bathtub: C = 0, n = 



Residues of Regge poles of draining bathtub: C = 0.5, n - 





FIG. 

(n 



8: Regge pole residues s^J^ in complex plane. The plots show the residues defined in Eq. 1 78 1 for the fundamental 
0) modes of a DBT, for the cases C — (non-rotating) [left plot] and C/D = 0.5 [right plot]. The blue and red 



crosses (solid lines) show the numerically-determined residues for 'prograde' (-I-) and 'retrograde' (-) modes, for frequencies 
LOVh = 1.0, 1.1, 1.2, . . . ,4.0. The dotted lines show the high-frequency approximation, Eq. 1 80 1. 



VI. FINAL REMARKS 

In this paper we have focused our attention on the 
quasinormal frequencies and Regge poles of the draining 
bathtub, a simple system which may be studied in the 
lab and which has been suggested as a simple (though 
inexact) analogue for the Kerr spacetime [17]. Our aim 
was to bring together a variety of techniques, both old 
[3T] and new p5H5Dl , to build up a unified picture 
of the resonances of a simple non-spherically-symmetric 
system, for the first time. 

We studied the QN resonances numerically in both 



the time domain (Sec. Ill A ) and the frequency domain 
(Sec. IIIB I. For the former we developed a purpose-built 



finite-difference code, and for the latter we applied the 
continued-fraction method (developed in Ref. PT]). We 
provided a clear geometric interpretation (Sec. Ill E I for 
our numerical results, in terms of the properties of the co- 
and counter- rotating photon/phonon orbits of the space- 



time (Sec. II B). The geometric interpretation neatly ex- 



plains the similarities and differences between the DBT 
and rotating black hole QN spectra (see Fig. [3| . Further- 
more, the geometric interpretation may be leveraged to 
obtain an expansion of QN frequencies in inverse powers 
of TO, as we showed in Sec. |III F[ 

In Sec. 



IV 



we demonstrated how methods developed 
to investigate QN modes may be simply adapted to inves- 
tigate the Regge pole spectrum (see also e.g. [15]). The 
relevance of the Regge poles is revealed by the Complex 
Angular Momentum method |47j . which we explored in 
Sec. |V] Regge poles may be used to account for oscil- 
lations seen in the absorption and scattering cross sec- 
tions. First, in Sec. V A we extended a method of [?7l[^ 



which reveals the role of Regge poles in fine-structure os- 
cillations in the absorption cross section (Tabs- We be- 
lieve this represents the first application of this work 
to a non-spherically symmetric spacetime. We showed 
(Fig. [7| that the superposition of diffraction effects linked 
to the properties of co- and counter-rotating orbits cre- 
ates a rather irregular fine-structure |28| in (Tabs('-^) when 
the system is rotating. Next, in Sec. |VB[ we applied 
the CAM method to obtain a geometric approximation 
for the 'orbiting' oscillations in the scattering cross sec- 
tion. To make this possible, we used the geodesic expan- 
sion method (56) to develop new approximations for the 
residues of the scattering matrix (Appendix |A]) , which 
are valid in the high-frequency limit. 

Let us bricfiy highlight some important analytic results 
herein: (i) at late times, perturbations of the DBT un- 
dergo power-law decay with an index 77 given by Eq. ( 20 1 ; 
(ii) the isotropic (to = 0) mode possesses a simple ana- 
lytic solution, Eq. (28 1 and (31); (iii) asymptotic expan- 



sions for QN and RP frequencies in terms of geomet- 
ric quantities are given by Eq. ( [5l| ) and (54), and they 
are good approximations in the large-|TO| and large-|a;| 
regimes, respectively (Fig. |4]); (iv) the absorption cross 
section may be expressed in terms of Regge poles via 
Eq. (64); (v) at large frequencies, the simple 'geometric' 



formula Eq. ( 71 ) provides a good approximation to Cabs 
(see Fig. [7]); (vT) the scattering amplitude at large angles 
is approximated by the Regge pole formula ( 77 ) ; (vii) in 



the high-frequency limit, the RP residues [required for 
(vi)] are given by Eq. (80). 



In conclusion, we have illustrated that a simple ana- 
logue system (the DBT) possesses a spectrum of reso- 
nances which, in addition to being of interest in their 
own right, can improve our understanding of the reso- 
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^ 0.12 - 




Scattering angle 



nances of a rotating black hole. Whereas black holes are 
remote celestial entities, analogues may be investigated 
in the lab today [TS] . We hope that a laboratory investi- 
gation of the resonances of the DBT will soon be under- 
taken. Here, experimentalists face the challenge of main- 
taining the stability of a converging flow, particularly at 
the point at which the flow becomes supersonic. It may 
be easier to instead use a diverging supersonic flow, which 
would instead represent a white hole analogue. For exam- 
ple, an hydraulic jump [16] extended to a rotating set-up 
may be of interest. It has to be noted, though, that 
the boundary conditions used to determine resonances 
of a white hole are different from the ones used in the 
corresponding black hole case [57l |58] . A further avenue 
for investigation is the effect of high-frequency dispersion 
on the resonance spectrum. For example, it has been 
argued that superluminal dispersion creates continuous 
zones in the QN mode spectrum of Bose-Einstein con- 
densates ^59j. Dispersion would certainly affect the va- 
lidity of the short-wavelength (i.e. large-m and large-w) 
approximations developed here, and would presumably 
disrupt the direct link between the properties of orbiting 
geodesies and the spectrum of resonances. Of further in- 
terest are 'pinned vortices' in superfluids, which resemble 
(to a certain extent) three-dimensional but non-draining 
versions of the DBT; it has recently been suggested that 
here bound states propagating along the vortex (similar 
to 'whispering gallery' modes) may arise |60| . 

To summarise, the key message from the present work 
is that the presence of distinct prograde and retrograde 
null geodesic orbits on a rotating (effective) spacetime 
leaves its clear imprint on the resonance spectrum. In 
turn, the superposition of co- and counter-rotating reso- 
nances generates a host of potentially-observable effects, 
for example (i) in QN ringing; (ii) in the fine-structure 
of the absorption cross section; and (iii) in the 'orbiting' 
diffraction effect in scattering. Here, we have applied a 
range of methods to investigate (i)-(iii) for the rather 
'idealized' geometry of the draining bathtub. It is our 
hope that these methods can now be extended to inves- 
tigate rotating analogue systems of current experimental 
interest. 



FIG. 9: Orbiting and the CAM approximation. These plots 
show the scattering cross section da/dcf) — \fuj{4i)\^ as a func- 
tion of scattering angle, with a coupling ujrh ~ 4, for three 
cases: C — (non-rotating), C = 0-2rh and C = l.Orh. 
The plots compare the cross section found from the partial 
wave series [red, solid] with the cross section found from the 
Regge pole (RP) approximation (77 1 [dotted line], using only 
the lowest modes n = 0...4). The blue dotted line shows 
the RP approximation with numerically-determined RPs and 
residues. The magenta dotted line shows the RP approxima- 
tion using asymptotic expressions ( 54 1 and 1 80 1 , for n = 
only. 



Appendix A: Regge pole residues 

In this section we apply the geodesic expansion method 
to obtain an asymptotic approximation for the Regge 
pole residues which is valid in the high-frequency limit. 
The method follows the outline given in Appendix A 
of Ref. [56] , where a similar calculation was performed 
for the (spherically-symmetric) Schwarzschild black hole. 
Here, the situation is complicated somewhat by the ro- 
tation of the system (C ^ 0). 

Recall that the residue is defined by Eq. (78). Here 
we will use the overbar notation fh = mj„ to represent 
the Regge pole value. To find the residue, we perturb 
the azimuthal number slightly away from its Regge pole 
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value, m = fh + e, (where e is small) to determine the 
first-order change in A™^. Away from the RP value, the 
'global' wavefunction ansatz (431 is no longer valid: its 
continuity breaks down near r « re- Instead, we may find 
regular solutions in a interior region close to r = Tc and 
match these onto exterior solutions (in regimes r < Vc 
and r > r^), to obtain = C'(e) to first order. In 

other words, we combine the geodesic expansion method 
with a standard WKB approach. 



Interior solution 



Let us begin with the radial equation (25) and first 
make the substitution u'^ = f~^/^X{r), to obtain 



where 



Uir) = /- 



Cm 



+ U{r)X = 0, 



(Al) 



1/4 7D'^ 



4r4 



^6 



Now, to find a solution valid in the vicinity of r 



make the change of variables r 



'^/"^z, where 



Substituting ( A2 1 into (All, with the assumption that lo 



is large, brings us to the following equation, 



^ + - 2^ (n + 7? + 1/2) + O ('^"'/')] = 0, 

(A3) 

where r\ = —itrdlc- Note that in our derivations we 
make use of the following identities: — CI 

2{rl - D^) = 2(1)2 _ (ji^y rpj^g solutions of Eq 

are the parabolic cylinder functions -D„+^((— l-|-i)z) and 
Dn+7j{{i~i)z). We use only the former solution, because 
(as we shall see) it matches onto an 'exterior' solution 
which is purely ingoing at the horizon. That is, we choose 



2 7' 



(A3) 



(A4) 



2. Exterior solutions 

In the exterior regimes r < Tc and r > rc, we require 
a pair of solutions (r) which are (i) approximate so- 
lutions to the radial equation in the large-w limit, and 
(ii) have well-defined 'ingoing'/'outgoing' behaviour in 
the limits — > ±cx). To obtain these solutions, we start 



with an ansatz of the form ( 43 ) , specifically, 



(r) exp 



a{r')f-\r')dr' 



(A5) 



with a as defined in Eq. ( |44| , and x^('') of the form 

exp(5±(r)). (A6) 



1 



To find explicit solutions, we expand the exponent a at 
orders (i) and (ii) and (iii) obtain by solving an 
equation like (49). After these steps (and finding some 



cancellation between terms at stages (ii) and (iii)), we 
reach the result 



±N-l/2 



r- D 
r + D 



( A2) where w = w - mC/D^ . 



(A7) 



3. Matching 

Let us now introduce three solutions, 

uq = Ai+»j((-l + i)^), r - Tc (A8) 
u< = r < r^. 

where and are coefficients to be determined. 

To match these solutions together in the transition re- 
gions, we make use of the following asymptotic forms: 



"""^ [2"/2e3""/4|^|«e+-V2+,^r(n+l)(2^)i/2|z|-(«+i)2-("+i)/2e-3'-("+i)/4e-»^'/2^ z ^ +oo. ^ ' 



and 



with A = 1/2, B = p/[2rca;i/2], and 



AF± {B\z\) 



±N-l/2 ±tz^/2 

^ 1 



(AlO) 



= exp ±i 



CjD ^ frc-D 

ujrc In 

2 \rr + D 



(All) 
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By matching the asymptotic forms, it fohows immedi- 
ately that 



A out 

^in 



a+ _ (^-)25-2W2We3.xJV/2 
^ r/(27r)i/2r(n+ 1) 



(A12) 



This leads us to the final (approximate) result for the 
residue, given in Eq. (80). 
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